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I 

INTRODUCTION 


The  stability  of  the  many  components  of  the  propblsion  system  must 
be  considered  jointly  In  order  to  establish  the  overall  stability  of 
the  system*  Because  of  the  nature  of  the  problem,  events  which  In 
themselves  take  a  minor  percentage  of  the  overall  stability'  margin,  cpn 
and  often  do  cause  stall t  Consequently,  it  is  all  too  easy  In  subse¬ 
quent  stall  Investigations  to  expend  more  effort  working  on  the  event 
which  tripped  the  stall  and  may  be  minor,  rather  than  the  more  funda¬ 
mental  phenomena  which  have  used  up  ninety-nine  pblnt  nine  percent  of 
the  stability  margin.  To  avoid  wasting  time  on  the  various  straws 
which  can  break  the  camel's  back,  it  is  necessary  to  develop  a  thorough 
accounting  system,  which  will  track  quantitatively  every  destabilizing 
parameter  from  the  early  stages  of  development  until  the  aircraft  is 
withdrawn  from  service. 

Computer  simulation  of  propulsion  components  and  systems  is  essential 
.  to  the  analysis  of  stability  margins.  The  math  models  for  these  simu¬ 
lations  must  however  rely  on  test  data  for  many  of  the  component  char¬ 
acteristics,  and  simulation  predictions  must  be  checked  according  to  a 
systematic  development  plan.  The  primary  role  of  the  computer  math  '\ 
model  is  to  synthesize  the  actions  and  interactions  of  the  many  compo¬ 
nent  parts  into  an  aggregate  system  response  for  a  specified  situation. 

The  aggregate  system  for  the  required ‘ stability  evaluation  may  comprise 
a  compressor,  an  engine  with  its  controls,  an  inlet  with  its  controls, 
the  complete  propulsion  system,, or  the  propulsion  system  with  a  flight 
dynamics  model  of  the  aircraft.  All  of  these  systems  levels  have  been 
modeled  for  a  variety  of  stability  analysis  purposes. 
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PROPULSION  SYSTEM  SIMULATION 


The  aoet  important  Models  for  compatibility  design  and  evaluation 
are  the  dynamic  Models  for  the  engine  with  controls  and  the  engine/ 
inlet  combi net Ion  (with  engine  and  Inlet  controls).  It  la  with  these 
Models  that  all  factors  contributing  to  surge  May  be  brought  together 
in  the  right  combination  and  order  dictated  by  the  imposed  disturbance, 
controls  limitations  and  tolerance  stacking  effects.  The  applications 
of  these  models  to  controls  design  and  node  selection  is  of  iaportance 
because  of  the  effect  which  the  controls  have  on  the  dynamic  mismatch 
between  engine  and  inlet.  Of  course  the  mismatch  is  important  to  both 
the  engine  and  inlet  stability  marginal  mismatch  is  often  a  source  of ^ 
distortion  and  turbulence  at  the  compressor  and  may  cause  buzz  or 
unatart  of  the  inlet. 

''  V  •  , 

Dynamic  simulation  of  engine /control  combinations  can  be  utilized 
to  define  operation  of  the  engine  and  to  assess  the  stability  margins 
and  controls  schedules  during  transient  excursions.  They  are  normally 
used  by  engine  manufacturers  during  a  development  program  to  develop 
the  controls,  monitor  component  requirements  and  to  coordinate  systems 
design  with  the  aircraft  manufacturer (a ) .  The  simulations  can  be 
developed  to  predict  propulsion  system  responses  to  all  of  the  required 
transients  (rotor  speed,  augmentation  changes  and  all. inlet  generated 
or  propagated  disturbances).  The  Inlet  generated  or  propagated  dis¬ 
turbances  (from  aircraft  maneuvers,  wind  gusts,  weapon  firing,  etc..) 
Include  steady  state  spatial  distortion,  time-variant  distortion,  inlet 
temperature  distortion,  and  Reynolds  number  effects.  To  evaluate  these 
latter  conditions  and  any  high-frequency  multi-dimensional  flow  distur¬ 
bances^  predetermined  correlations  in  component  performance  and  stabil¬ 
ity  characteristics  must  be  included  in  the  mathematical  model.  Test 
data  and  other  analytical  procedures  can  be  used  to  obtain  the  required 
correlations  between  distortion,  etc.  and  the  changes  in  component 
performance  and  stability  margin. 

/ 

The  use  of  a  nonlinear  dynamic  model  of  the  engine  and  control  has 
gained  wide  acceptance  In  the  industry.  This  technique  requires  the 
definition  of  the  engine  component  performance  and  stability  character¬ 
istics  in  the  form  of  compressor  maps,  turbine  maps,  burner  stability 
limits,  etc.  The  performance  maps  of  each  of  the  components  can  be 
modified  to  reflect  the  effects  of  distortion  and  turbulence.  Figure  1 
describes  the  performance  information  required  for  a  split  compressor 
turbo fan  engine.  Each  of  the  components  is  treated  as  a  lumped  element, 
and' major  gaa  storage  elements  are  treated  as  lumped  volumes.  The  choice 
of  volume  lumpings  shown  by  Figure  2  is  for  illustrative  purposes  only, 
the  final  number  of  volumes  chosen  is  a  function  of  the  frequency 
response  required. 
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The  dynamic  coupling  beCveen  components  in  the  engine  simulator  is 
accomplished  by  two  basic  relationships:  (1)  rotor  speed  of  each  spool 
is  obtained  from  the  tlM  Integration  of  the  imbalance  In  torque 
(divided  by  rotor  inert la)  b^tHsen  the  turbine  and  compressor  on  that 
rotor;  and' (2)  pressure  between  components  is  obtained  by  time  Integra-^ 
tion  of  the  airflow. Imbalance  (as  modified  by  volume,  specific  heat 
ratio,  and  temperature  according  to  the  ideal  gas  law)  between  inter¬ 
faces  of  the  control  volume.  The  more  precise  models  for  higher  fre¬ 
quency  must  incorporate  energy  conservation  equations  (especially  for 
burner  and  augmentor)  involving  enthalpy  flows  at  interfaces  and  dynamic 
combustion  equations.  The  computerised  dynamic  math  model  indicated  in 
Figure  1,  with  variations  ss  described  above,  la  useful  in  the  study  of 
stability  margin  variations  resulting  from  the  interaction  of  provisional 
controls  definitions  and  s'  variety  of  system  disturbances.  The  dynamic 
reeponae  for  this  type  of  model,  can  usually  be  made  apod  to  approximately^ 
20, hertz.  Effect^  requiring  higher  frequency  response  must  be  accounted 
fdr  by  modification  of  the  basic  amps. 


/ 


/ 


/ 

/ 

/ 

/ 

/ 
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III 

STABILITY  PREDICTION  VALIDITY 


Propulsion  system  simulations  have  demonstrated  the  capability  to 
follow  actual  engine  transients*  The  mathematical  capability  and  under¬ 
standing  of  the  relevant  thermodynamic  relationships  can  therefore  be 
considered  to  be  in  hand.  However,  the  margin  as  assessed  by  the  simu¬ 
lation  will  only  be  accurate  provided  the  input  empirical  data  is 
correct.  Some  of  the  data  input  requirements  are  straightforward,  for 
example,  component  performance  maps,  engine  volumes  and  areas,  control 
schedules  and  component  Inertias.  1“ fortunately,  some  of  the  required 
data  is  less  straightforward;  thiamin  general  consists  of  the  modifi-  1 
cation  necessary  to  the  basic  maps  forced  by  component  interactions 
and  the  effect  of  changes  in  external  environment.  It  is  the  degree 
to  which  these  interactions  and  modifications  can  be  accounted  for  that 
nukes  the  difference  between  a  simulation  which  'is  an  Interesting  toy 
and  a  simulation  which  is  a  powerful  and  practical  technique  for  keeping 
track  of  propulsion  system  stability. 

-Past  experience  can  be  used  initially  for  factors  such  as  Reynolds 
Number  effects  on  compressor  performance,  but  this  can  be  such  a  signi¬ 
ficant  effect  that  subsequent  test  verification  is  required.  Some  map 
changes,  however,  may  require  a  major  test  program  to  support  the 
necessary  modifications.  Figure  3  shows  s  basic  fan  map  as  tested  on  a 
rig.  Figures  4  and  S  show  the  changes  In  the  map  of  the  same  fan  when 
the  splitter  cowl  which  separates  engine  core  flow  from  the  fan  duct 
flow  is  placed  close  to  the  fan  and  the  two  streams  are  separately  back- 
pressured.  Figure  4  shows  the  change  ,in  the  map  when  the  fan  duct 
stream  is  back-pressured.  Figure  5  shows  the'  change  In  the  fan  sup 
when  the  core  stream  la  back-pressured.  In -actual  engine  operation,  not 
only  will  the  surge  line  change  as  a  function  of  engine  operating  con¬ 
dition,  but  the  operating  line  will  he  similarly  effected.  Clearly  a 
modification  of  this  type  to  the  simulation  is  a  difficult  programsilng 
job.  However,  in  the  long  run  it  (is  preferable  and  less  prone  to  error 
than  attempting  to  track  by  hand  the  continuously  changing  margin 
through  an  engine  transient ,  particularly  since  all  other  destabilizing 
Influences  must  be  monitored  at  the  same  time. 

i 

J 

For  many  applications  the  key  parameter  in  determining  propulsion 
system  stability  is  the  index  used  to  describe  the  flow  field  at  the 
engine  face.  The  history  of  distortion  Indices  has  been  somewhat 
checkered.  Many  early  indices  were  claimed  by  their  inventors  to  be 
applicable  to  all  engines.  This  was  a  fallacy  but  index  Inventors 
often  stuck  by  their  numbers  as  though  their  personal  integrity  and 
their  company's  position  depended  on  proving  the  index  vas  valid.  i 
Figure  6  and  Figure  7  show  comparisons  of  two  pairs  of  distortion 
indices.  The  pair  in  Figure  6  includes  a  dynamic  head  (Q)  term  and  the 
pair  in  Figure  7  do  not.  Each  point  represents  a  different  distortion 
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pattern.  An  index  of  a  particular  value  ia  supposed  to  represent  a 
given  increaent  of  aurge  preaaure  ratio  lose.  It  ia  clear,  therefore, 
that  if  the  indicea  are  applicable  to  ell  engines  a  relationship  should 
exist  between  the  indices  which  should  be  single  valued.  This  single 
valued  relationship  does  not  exist  unless  a  given  form  of  pattern  is 
specified,  hence  one  index  la  not  generally  applicable  and  the  probabil¬ 
ity  is  that  neither  are.  The  situation  was  further  confused  by  the  fact 
that  in  the  early  days  distortion  on  aircraft  waa  measured  using  steady 
state  instrumentation.  Hence,  even  when  an  Index  was  tailored  to  a 
particular  engine,  there  was  no  correlation  with  flight  test  data 
because  the  engine  was  seeing  distortion  patterns  which  were  quite 
different  from  those  recorded  by  the  instrumentation.  At  this  point 
the  index  inventors  were  accused  of  running  a  numbers  game.  This  feeling 
is  in  no  way  dissipated  when  some  of  the  "confusion  factors"  which  go 
into  indices  are  investigated.  As  an  example,  consider  the  use  of  the 
dynamic  head  Q.  This  is  claimed  to  be  a  normalising  factor,  since  it 
compensates  at  part  speed  for  the  fact  that  a  given  level  of  delta 
pressure  corresponds  to  a  greater  change  in  velocity  than  that  same  delta 
pressure  at  high  speed. 

Figure  8  shows  the  loss  of  surge  pressure  ratio  of  the  fan  as  a 
function  of  speed  (i.e. ,  flow)  and  a  distortion  index.  Figure  9  shows 
the  same  data  except  that  this  time  the  distortion  index  has  been  norma¬ 
lized  (l.e. ,  allegedly  compensated  for  change  in  flow)  by  use  of  a  Q/Pt 
factor.  More  interesting  games  can  be  played  using  this  factor  when  con¬ 
sidering  fan  amplification  or  attenuation  of  distortion.  Figure  10  shows 
a  distortion  transfer  coefficient  (distortion  out/dlstortlon  in)  plotted 
in  terms  of  an  index  based  on  delta  pressure  in/delta  pressure  out  and 
the  same  index  normalised  by  a  Q  factor.  In  this  cash,  two  different 
stations  were  used  for  calculating  the  input  Q.  Although  the  component 
performance  is  the  same  for  each  plot,  the  apparent  performance  to  the 
uncritical  eye  can  be  greatly  improved  by  aui table  choice  of  index. 

Deapite  the  prohleaui  aaaoclated  with  dlatortion  indicea,  it  la  clear 
that  the  flow  field  at  the  compressor  face  must  be  described  by  a 
numerical  syatem  if  quantitative  assessment  of  available  margin  is  to  he 
achieved.  Bearing  in  mind  the  apparent  component  improvement  implied  by 
the  previous  example,  it  is  extremely  apparent  that  any  system  of  sta¬ 
bility  resolution  is  dependent  on  complete  visualization  of  the  flow 
processes  of  each  of  the  components  Involved.  However,  organization  of 
the  program  requires  the  Interfacing  of  many  technical  parameters  to 
form  the  quantitative  basis  for  assessment  of  stability  margin.  It  is 
essential,  however,  to  resiember  and  watch  for  the  pitfalls  associated  with 
distortion  Indices.  The  man  who  offers  to  sell  distortion  indices  and 
the  man  who  offers  to  sell  the  Brooklyn  Bridge  stay  not  be  out  of  the  same 
stable,  but  it  would  be  as  well  to  check. 

The  first  step  in  establishing  a  workable  quantitative  relationship 
between  distortion  and  surge  margin  loss  is  to  establish  the  fan  response 
to  steady  state  distortion.  Several  functions  must  be  considered.  First, 
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sensitivity  ■—  the  rate  of  loss  of  surge  pressure  ratio  with  both  circum- 
ferentlal  and  radial  distortions  sust  be  found  from  rig  tests.  Distortion 
sensitivity  will  In  vary  float  hub  to  tip.  In  sons  cases,  this 

is  e  very  pronounced  effect  end  the  lnder  should  account  for  It.  The 
extent  and  nunber  of  lobes  of  a  circumferential  distortion  has  e  major 
effect,  also  the  prassure  gradient  as  the  blades  enter  and  leave  a  low 
pressure  sons  nay  be  significant.  The  loss  of  surge  pressure  ratio 
associated  with  these  phenomena  should  be  Investigated  using  controlled 
classical  patterns.  All  the  effects  should  then  be  added  together  end 
the  Index  system  checked,  using  more  complex  patterns  which  are  typical 
of  actual  Inlet  distortion.  At  the  some  time,  the  fan  transfer  coeffi¬ 
cients  (l.e.,  the  relationship  betwean  distortion  in  and  distortion  out) 
should  be  established.  Unless  the  distortion  into  the  compression 
system  behind  the  fan  can  be  accurately  determined  knowing  only  the 
conditions  in  front  of  ths  fan,  diagnosis  or  assessmsnt  of  overall  engine 
stability  is  extremely  difficult.  Furthermore,  It  is  impossible  to 
screen  Inlet  data  for  distortion  levels  potentially  dangerous  to  the 
engine  without  adequate  knowledge  of  the  fan  transfer  coefficients. 

These  coefficients  should  account  for  both  redial  distortion  and  the 
extant  depth  and  lobe  pattern  of  circumferential  distortion.  Zn  addition 
to  the  pressure  stability  coefficients,  It  Is  necessary  to  account  for 
the  temperature  distortion  which  is  introduced  by  the  different  fan  work 
inputs  in  the  high  and  low  pressure  zones  of  the  inlet  flow  field. 

The  stability  coefficients  for  the  subsequent  units  In  the  compres¬ 
sion  system  should  be  established  In  the  seam  way  as  for  the  fan.  However, 
two  additional  points  must  be  bourne  in  mind.  The  dean  surge  line  must 
be  established  with  the  radial  profile  Imposed  which  is  normally  delivered 
by  the  forward  unit  (or  units),  and  the  Imposed  distortion  must  be  cal¬ 
culated  after  applying  the  tranafer  coefficients  for  the  forward  units 
to  the  relevant  section  of  the  inlet  flow. 

Once  a  credible  numerical  description  has  been  developed  for  the 
inlet  flow  field,  it  Is  possible  to  proceed  with  snalysis  of  a  real 
dynamic  environment,  The  patterns  which  the  engine  is  encountering  or 
can  be  expected  to  encounter  must  ba  determined.  In  a  full  scale  test, 
the  relevant  patterns  may  be  found  by  examining  a  moving  time-average  of 
the  instantaneous  distortion  measured  by  high  response  instrumentation. 

The  test  date  end  instrumentation  requirements  on  which  this  analysis 
technique  is  based  are  described  in  Reference  1.  This  analysis  technique 
eliminates  the  error  which  vas  introduced  when  inlet-distortion  data 
was  based  on  the  readings  of  steady  state  or  very  low  response  instru¬ 
mentation,  but  there  are  other  problems  of  environmental  definition. 

When  only  inlet  model  data  is  available,  l.e.,  in  the  early  stages  of 
system  development,  great  care  must  be  taken  to  account  for  scale  effects. 
The  effects  of  model  size,  acoustic  frequency  scaling  and  Reynolds 
Number  must  be  applied.  A  fuller  discussion  of  dynamic  distortion 
scaling  may  be  found  in  other  references. 


6 


AFAPL-TR-70-57 


. 80  far  *n  tha  di.cuj.ion  of  th.  data  requirements  for  a  .tabilitv 
atack  up,  only  passing  rsfsrsncss  havs  bean  asde  as  to  the  tine  ohasiL 

2Le^.ST '  *2£T!S»;.  7“ *•  «w “f 


- - Will 

research  projections, 
■uet  be  replaced  by  hard 
whenever  subsequent  tests 


UB  Da*®°  011  tn*  «nglne  aanuf acturers '  experience 
As  the  program  develops,  these  projections 
test  rig  data  which  must  then  be  updated 
are  run. 


AFAPL-TR-70-57 


IV 

COBCUJSIOH 


It  is  currently  possible  to  combine  in  a  propulsion  system  dynamic 
model  the  empirical  and  theoretical  data  nacessary  to  determine  the 
available  transient  stability  margin  at  various  fiight  conditions^  The 
date  requirements  for  such  a  model  are  massive,  but  they  are  essential 
for  satisfactory  stability  development,  and  none  of  the  data  r«*lu^***“ 
necessitate  techniques  outside  the  current  state-of-the-art.  In  *“***“ 
to  stability  tracking,  the  model  provides  at  any  stage  in  the  Jeveiopmant 
program  a  guide  as  to  the  work  remaining  to  be  done  and  a  check  on  the 
progress  of  the  work  already  accomplished. 
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1  PERFORMANCE  MAPS  REQUIRED 
DYNAMIC  ENGINE  SIMULATION 


FIGURE  2  TYPICAL  LUMPED-VOLUME  ELEMENTS  AND  COUIPLMG 

FOR  A  DYNAMIC  ENGME  SIMULATION 
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FIGURE  3  TYPICAL  FAN  MAP 


FIGURE  4  EFFECT  OF  O.D.  BACK  PRESSURE 
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FIGURE  6  OVERALL  DISTORTION  INDEX  COMPARISON 

INCLUDING  Q  FACTOR 


■  180°  PATTERN 
•  COMPLEX  PATTERN 


CIRCUMFERENTIAL  DISTORTION  INDEX  COMPARISON 

NO  Q  FACTOR 


MARGIN  LOSS  AT  100%  SPEED  =  2 X  MARGIN  LOSS  AT  50%  SPEED 
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AND  SPEED  NORMALIZED  BY  Q 


Hill 


PERCENT  FLOW 

Q  ON  DISTORTION  TRANSFER  COEFFICIENT 


^SS^mS^ISSSSmSStSZJ^lSSSSSSiL 

RCPrRT  IICOHITV  CLASSIFICATION 


xnajkSSiruD 

hcuili  CUiciAcctlM 


STABILITY  MARGIN  DEFIHITIOS 
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II.  IUMLIMIIITAIIV  NOTH 


IP.  APITRACT 


The  aval lab la  stability  nergln  of  a  turbine  propulsion  eye ten  la  a  function  of  a 
wide  range  of  parameters ,  natty  of  which  are  tine  dependant.  Both  Internal  and 
external  effects  are  considered  and  the  interactions  between  Inlet  and  engine 
stability  coefficients  are  discussed.  (  i 
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UNCLASSIFIED 
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STABILITY  MARGIN 


OVERALL  STABILITY  KARGIN 
ENG IKE /INLET  STABILITY 
STABILITY  COEmCIENTS 
TRANSIENT  STABILITY  MARGIN 
INLET  STABILITY 
TURBINE  ENGINE  STABILITY 
INLET  DISTORTION 
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